Adatom Diffusion at GaN (0001) and (0001) Surfaces 
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The diffusion of Ga and N adatoms has been studied for the technologically relevant wurtzite 
(0001) and (0001) surfaces employing density-functional theory. Our calculations reveal a very 
different diffusivity for Ga and N adatoms on the equilibrium surfaces: While Ga is very mobile at 
typical growth temperatures, the diffusion of N is by orders of magnitudes slower. These results 
give a very detailed insight of how and under which growth conditions N adatoms can be stabilized 
and efficiently incorporated at the surface. We further find that the presence of excess N strongly 
increases the Ga diffusion barrier and discuss the consequences for the growth of GaN. 



Recently, great progress in fabricating highly efficient 
GaN-based devices has been achieved Nevertheless, 
there are still substantial problems concerning growth op- 
timization and insight into the fundamental mechanisms 
is rather shallow. Investigations of molecular beam epi- 
taxy (MBE) growth have shown that the film structure 
and morphology are very sensitive to the III/V ratio ||f| . 
GaN films grown under more N-rich conditions are rough 
and faceted while going towards more Ga-rich conditions 
a smoother surface morphology and better film quality 
are obtained || . 

In order to improve growth in a systematic way it is 
essential to understand the underlying kinetic processes 
such as adsorption, desorption, and surface diffusion. In 
particular, adatom diffusion on surfaces is considered to 
be a key parameter controlling the growth rate, the ma- 
terial quality, and the surface morphology. Experimen- 
tally, an analysis of surface diffusion is difficult: So far 
only effective diffusion barriers have been obtained for 
GaN [0||. It is also not clear whether the cation or an- 
ion surface diffusion is the rate limiting process Q. A 
problem in growing GaN is also the high N vapor pres- 
sure and considerable efforts have been made to enhance 
the N incorporation. However, the mechanism by which 
N is incorporated at the surface could not be identified. 
Theoretically, only few studies to describe GaN growth 
have been reported, which were either based on ther- 
modynamic models || or Monte Carlo simulations 0. 
These approaches are useful to model growth on a meso- 
scopic scale. However, the microscopic meaning of the 
effective parameters and the underlying microscopic pro- 
cesses remain unclear. In order to identify the fundamen- 
tal growth mechanisms, but also to improve the above 
mentioned methods, the correct microscopic parameters 
are needed [§. 

We have therefore performed a comprehensive and 
detailed study of the migration and energetics of Ga 
and N adatoms on GaN surfaces, employing total- 
energy density-functional theory calculations with ab ini- 
tio pseudopotentials. We will focus on the two technolog- 



ically relevant orientations for wurtzite GaN, the (0001) 
and the (0001) surfaces. The adatom substrate system is 
modeled by supercell geometries with at least 9 layers of 
GaN, a 14 bohr vacuum region and (2 x 2) periodicity. 
One side of the slab is passivated by fractional pseudo- 
hydrogen. Tests revealed that migration paths and dif- 
fusion barriers are not significantly affected by the Ga 
3d electrons and that it is sufficient to treat them within 
the non-linear core correction (NLCCJjIJ. Details of the 
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calculations are discussed elsewhere 

The energetically stable reconstructions for polar GaN 
surfaces have been recently identified for the cubic (001) 
@ and the wurtzite (0001) and (0001) orientations 
These studies revealed a surprising feature not observed 
for "traditional" III-V semiconductor surfaces: Indepen- 
dent of the chemical environment (Ga- or N-rich) polar 
GaN surfaces turn out to be always cation stabilized with 
no N atoms in or on-top of the surface layer. The only 
known exception so far is a N (2 x 2) adatom structure 
on the (0001) surface found to be stable under very N- 
rich conditions Jl3[ ] . In order to study adatom diffusion 
on the thermodynamical stable surfaces we will therefore 
focus on the two structures shown in Fig. 1 which both 
are characterized by a complete Ga surface layer. 

We first calculated the potential energy surface (PES) 
for Ga and N adatoms on the (0001) surface, giving 
immediate insight into stable sites, migration paths, 
and diffusion barriers. The potential energy surface 
i?tot(R| d ) is calculated by fixing the adatom laterally at 
different positions R| d and allowing all other atoms and 
the adsorbate height to relax. A two-dimensional (2D) 
cut through the resulting potential energy surfaces for 
Ga and N adatoms along the energetically lowest path is 
shown in Fig. 2. For Ga adatoms we find two adsorption 
sites: fee and hep (see Fig. 1), which are degenerate in 
energy within the accuracy of our method. Fig. 2 also 
reveals that the energetically lowest transition site [corre- 
sponding barrier: (0.4 eV)] is the bridge position; and the 
migration over the on-top position is energetically very 
unfavorable. Actually, an analysis of the complete PES 
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showed that the on-top position is the global maximum. 

For N adatoms a qualitatively very different PES re- 
sults (Fig. 2, dashed line): While the fee position is 
the energetically lowest binding site, the hep position is 
significantly higher in energy (1.4 eV). Therefore, to hop 
from one fee site to the next, the N adatom diffuses along 
the bridge position over the hep site which is the transi- 
tion site leading to a barrier of 1.4 eV. The large energy 
difference between the fee and hep sites can be under- 
stood by analyzing the atomic structure of the surface: 
At the hep position, the N adatom experiences directly 
underneath a N atom in the second layer. Since GaN is a 
partly ionic material this leads to electrostatic repulsion 

B 

When analyzing adsorption sites it should be kept in 
mind that the on-top position is the epitaxial site for N 
adatoms on the (0001) surface, i.e., for high coverages this 
site must be a global minimum. From Fig. 2 we find this 
site to be energetically highly unfavorable for an isolated 
N adatom (i.e. for low coverages). Therefore, when going 
from low to high adatom coverages the N adatoms have 
to shift from the fee site (stable at low coverages) to the 
on-top position (stable at high coverages). Since this 
process is kinetically hindered, the initial adsorption at 
the fee position might be a potential source for stacking 
faults. 

We have also studied the diffusion on the (0001) sur- 
face, which is characterized by a very different atomic 
geometry where each surface atom has three dangling 
bonds compared to only one at the (0001) surface (see 
Fig. 1). For Ga adatoms the PES (see Fig. 3) shows 
very similar features as found on the (0001) surface: The 
fee and hep sites are energetically degenerate, the diffu- 
sion transition state is at the bridge site with an again 
very low energy barrier (0.2 eV). 

For N adatoms we find a qualitative different behav- 
ior compared to N atoms on the (0001) surface (see Fig. 
3): (i) N adatoms exhibit the same binding sites than 
Ga adatoms and (ii) the diffusion barrier is significantly 
lower [0.9 eV compared to 1.4 eV on (0001)]. The lower 
diffusion barrier can be understood in terms of the very 
different character of the chemical bonds on the (0001) 
and (0001) surfaces. On the (0001) surface (see Fig. 1) 
the localized and strongly directed sp z orbitals favor the 
formation of strong Ga-N bonds giving low coordinated 
sites, as e.g. the bridge position, a rather high energy. 
On the other hand, the (0001) surface is characterized by 
metallic bonds between the surface atoms: the adsorbate- 
substrate interaction is hence significantly weaker result- 
ing in lower diffusion barriers. 

An important consequence of our results is that Ga 
adatoms have a significantly lower diffusion barrier than 
N adatoms for both orientations: Ga adatoms will be 
orders of magnitude more mobile than N adatoms at 
typical growth temperatures. The low diffusion bar- 
rier is a direct consequence of the fact, that for GaN 



equilibrium surfaces are Ga stabilized. Thus, for Ga 
atoms the adsorbate-substrate interaction is predomi- 
nantly realized by delocalized metallic Ga-Ga bonds. 
Since Ga-bulk melts already slightly above room temper- 
ature (Tmcit = 30°C), the Ga-Ga bonds are weak and the 
adatoms behave almost like a liquid film on the surface. 
A similar effect has not been reported at "traditional" 
III/V semiconductor surfaces like e.g. GaAs where the 
surfaces do not exhibit a metallic like character. The 
diffusion barrier on these surfaces is thus mainly char- 
acterized by breaking strong cation-anion bonds. As a 
consequence a significantly higher Ga diffusion barrier is 
found on these surfaces, e.g., for the polar GaAs (111) 



surface the barrier is 0.9 eV 15|. 

The significantly larger diffusion barrier for N adatoms 
compared to Ga adatoms on equilibrium surfaces has an 
important consequence: Although N adatoms on these 
surfaces are thermodynamically unstable against evapo- 
ration as N2-molecules even under N-rich conditions jl2) , 
they can be kinetically stabilized at the surface. In order 
to evaporate, two N atoms have to form a N2 molecule. 
Since migration of N adatoms is a highly activated pro- 
cess (but necessary to form molecules), the desorption 
rate may become smaller than the adsorption rate. Con- 
sequently, extended regions in which the surface is pri- 
marily covered by N atoms may be formed. These N 
adatoms are likely to influence the migration path and 
the diffusion barrier of Ga adatoms. We have therefore 
studied the diffusion of Ga adatoms on the N-terminated 
(0001) and (0001) surfaces which can be considered as 
"extreme" cases of N coverage. The results are shown 
in Fig 4. The energetically favored binding sites are lo- 
cated at the fee and hep positions and the transition site 
is the bridge position. For both surface orientations the 
diffusion barrier of Ga adatoms is strongly affected: At 
(0001) the migration barrier increases from 0.4 to 1.8 eV 
while at (0001) it increases from 0.2 to 1.0 eV. We there- 
fore conclude that excess N at the surface, which can be 
formed under N-rich conditions, significantly reduces the 
mobility of Ga adatoms. The reason is the formation of 
strong Ga-N bonds, which have to be broken during the 
migration of the adatoms. 

The reduction in surface diffusivity when going to- 
wards more N-rich conditions is consistent with a recent 
very detailed MBE growth study by Tarsa et al. S. We 
note however, that a direct comparison between our cal- 
culated diffusion barriers and experimentally derived bar- 
riers is by no means straightforward. Experimentally, ef- 
fective diffusion parameters are commonly derived which 
average over a large area including not only the clean sur- 
face but also steps, impurities, dislocation etc. Since the 
effective diffusion barrier is usually dominated by the pro- 
cess with the highest barrier, our calculated value (which 
is for the clean surface) is a lower limit of the effective 
diffusion barrier. 

The very different mobilities of Ga and N adatoms 



2 



and the formation of excess N have important conse- 
quences for the growth of GaN. In the Ga-rich regime 
where the amount of excess N on the surface is small, 
the Ga adatoms are highly mobile and a step-flow mode 
resulting in 2D growth is expected. As a consequence, 
the surface morphology should be improved and a low 
density of stacking faults is expected. Further, if ex- 
cess Ga adatoms are present on the surface (as expected 
for Ga-rich growth conditions), N adatoms can be effi- 
ciently incorporated: The probability that fast moving 
Ga adatoms capture N atoms is much higher than the 
other process where N atoms form molecules and desorb 
from the surface. However, at the N-terminated surfaces 
our results show a roughly five times higher diffusion 
barrier, indicating that the Ga diffusion length is sig- 
nificantly shorter under N-rich conditions. Once the dif- 
fusion length is shorter than the mean distance between 
the binding sites a statistical roughening of the surface 
can be expected. Furthermore, the adatoms might be 
"trapped" at sites not corresponding to the ideal bulk 
positions. For example, if a Ga atom is "trapped" at 
the fee position and does not hop to the wurtzite hep 
site, a fee nucleation center is formed. Thus, the higher 
adatom mobility under more Ga-rich conditions will also 
significantly reduce the density of stacking faults. Based 
on a detailed analysis of the adatom kinetics we expect 
therefore slightly Ga-rich conditions to be optimal for the 
growth of GaN. We note however, that while going to- 
wards extreme Ga-rich conditions might further improve 
the growth morphology, it will also enhance the impurity 
incorporation || and result in high background carrier 
concentrations. 
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FIG. 1. Atomic structures (top and side view) for the 
(0001) and (0001) Ga terminated surfaces. The dashed line 
in the top left marks the diffusion pathway corresponding to 
Fig. 2 and Fig. 3. Ga±L marks a Ga-atom in the first layer, 
Ga2L in the second. 
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FIG. 2. Total energy (in eV) for Ga-adatom (solid line) 
and a N-adatom (dashed line) at the Ga-terminated (0001) 
surface. The energy zero corresponds to the energetically low- 
est adsorption site. 
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FIG. 3. Total energy (in eV) for a Ga-adatom and a 
N-adatom at the Ga-terminated (0001) surface. The energy 
zero corresponds to the energetically lowest adsorption site. 
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FIG. 4. Total energy (in eV) for Ga-adatoms at the 
N-terminated (0001, solid line) and the N-terminated (0001, 
dashed line) surface. The energy zero corresponds to the en- 
ergetically lowest adsorption site. 
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